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Abstract 

We prove that for divergent series solutions of nonlinear (or linear) 
differential systems near a generic irregular singularity, the common pre- 
scription of summation to the least term is, if properly interpreted, mean- 
ingful and correct, and we extend this method to transseries solutions. In 
every direction in the complex plane at the singularity (Stokes directions 
not excepted) there exists a nonempty set of solutions whose difference 
from the "optimally" (i.e., near the least term) truncated asymptotic se- 
ries is of the same (exponentially small) order of magnitude as the least 
term of the series. 

There is a family of generalized Borel summation formulas B which 
commute with the usual algebraic and analytic operations (addition, mul- 
tiplication, differentiation, etc). We show that there is exactly one of 
them, Bo, such that for any formal series solution /, Bo(f) differs from 
the optimal truncation of / by at most the order of the least term of /. 

We show in addition that the Berry (1989) smoothing phenomenon 
is universal within this class of differential systems. Whenever the terms 
"beyond all orders" change in crossing a Stokes line, these terms vary 
smoothly on the Berry scale arg(a;) ~ |a;|~ 1/ ' 2 and the transition is always 
given by the error function; under the same conditions we show that 
Dingle's rule of signs for Stokes transitions holds. 

1 Introduction 

Summation to the least term (optimal truncation of series) can be traced back 
to Cauchy's study of the Gamma function (see Cauchy (1882)). The method 
was applied by Stokes (see Stokes (1904)) to solutions of differential equations 
and was instrumental in his theory of Stokes' phenomenon. Optimal truncation 
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techniques have been greatly improved by Dingle (1973), and in recent years 
as a result ol the new ideas and methods ol hyperasymptotics introduced by 
Berry (1989), who also discovered a surprisingly universal transition in hyper- 
asymptotics, often called Berry smoothing . Rigorous results on optimal trunca- 
tion and smoothing under various assumptions were proved by McLeod (1992) 
(second order linear equations), Paris (1992) (linear systems with power coef- 
ficients), Olvcr (1995) (Gamma function), Olver and Olde Daalhuis (1995a, b) 
(linear second order systems) , Olde Daalhuis (1995, 1996) (linear second or- 
der systems), the authors (1996) (first order nonlinear and second order linear) 
and others. A wide range of interesting applications were considered by Berry 
and Howls (1990, 1993) (saddle point integrals), Berry and Keating (1992) 
(asymptotics of the zeta function) and Berry (1991, 1994, 1995) (asymptotics 
of the Gamma function, applications to dynamical systems, and PDE's). Other 
techniques of exponentially improved asymptotics have found interesting appli- 
cations in dynamical systems (see Neishtadt 1984, Ramis and Schafke 1996). 

In this paper we prove that asymptotic series of solutions of generic ana- 
lytic (linear or nonlinear) differential systems near irregular singularities are 
summable to the least term, with errors of the order of magnitude of the least 
term with respect to the functions associated to the series by a specific gener- 
alized Borel summation (the balanced average). In fact we show that the same 
is true for all series components of the transseries solutions of such differential 
systems. 

We show that Berry smoothing extends to this class of asymptotic series. 
We deal with equations of the form 



y' = f(*,y) yeC" (1.1) 

for large x in some direction, and study solutions that go to zero in this limit. 
The assumptions on H.lfl arc 

(al) The function f is analytic at (oo,0). 

(a2) The eigenvalues Xi of the linearization 

k:=-(&oo,0)) (1.2) 



9yj 



l,2,...n 



arc different from zero and nonresonnant (cf. ijl.lf) . 

It is worth mentioning that a large class of differential systems are not pre- 
sented in the form described above but can be brought to it by suitable changes 
of variables. For example, under the change of variables 2a; 3 / 2 /3 = t, the Airy 
equation y" — xy = becomes 

d 2 y , 1 dy 

Also, the Painleve I equation y" = 6y 2 + x becomes 
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d 2 h ldh , 3,, 392 

— H h - -h 2 r = 

dt 2 t dt 2 1875i 4 

by taking y(x) = -i\/6i/2[l/3 - 4^ 2 /75 + h(t)} and t = (-24x) 5 / 4 /30. In 
this form, written as systems, both equations satisfy our assumptions (see HI. II 
and the notes therein). The transformations themselves, when they exist, are 
readily obtained by comparing the complete formal solution (transscrics) of the 
original equation with the transseries of the normalized equation, (|1.4|) below. 

It is convenient to pull out the inhomogencous and the linear terms (relevant 
to leading order asymptotics) and rewrite the system in the form 

y' = f (x) - Ay - -By + g(x,y) (1.3) 

x 

Under the assumptions (al) and (a2) equation (|1.3|l admits, in normalized 
form, an n— parameter family of formal exponential series solutions 



y = yo+ Yl cJ ? •••^"^ (kA) ^ km y k (1.4) 

k>0;|k|>0 

(see HI .2(1 where yk = x~^P +ra ^ X^o ak ; iX_ ' are f° rma l power series. 

We study those formal solutions (|1.4|l which are at the same time asymptotic 
expansions (or proper transscrics, in the sense of Ecalle 1993). In our context, 
only finitely many powers can be present, of any exponential which is not small 
in the given direction. Note however that any exponential will decrease in some 
direction, and our analysis eventually covers an n parameter family of formal 
solutions. 

For linear equations, any formal solutions l|1.4f) contain finitely many expo- 
nentials and are proper transseries. For nonlinear equations we require, for large 
x in a given direction d in C, that Cj = for any i such that e~ XiX -/-* (more 
details in HT^ . 

We start by illustrating the various concepts on a very simple example. 
The equation /' + / = 1/x has the general solution y{x) = e~ x Ei(x) + Ce~ x 
(Ei(a;) := PV J Q t~ 1 e t dt). The general formal solution for large x is the el- 
ementary transseries yc = SfcLo klx~ k ~ 1 + Cc~ x with C G C arbitrary. If 
arg(x) € (7r/2, 37r/2), the exponential term is large and classical asymptotics 
gives full meaning to the family yc'- for each yc there is a unique actual solu- 
tion yc of the equation such that yc ~ yc f° r large x. But as arg(x) crosses an 
antistokes line arg(x) = (fc + l/2)7r into a sector where e~ x decays, e~ x becomes 
small beyond all orders of the divergent series and is undefined (as part of yc) 
in classical asymptotics. Classical asymptotics gives up the exponential and re- 
tains the information that all solutions are asymptotic to yo. In this framework 
there is no natural way for asymptotically choosing any privileged solution as 
corresponding to yo; yo becomes more of a property of the differential equation 
as a whole, shared by all its particular solutions. 



3 



In terms of estimating solutions out of yo, classical asymptotics provides 
polynomial precision: \y(x) — < Const. la;! - ^" 2 for large x. where y N ^ := 

One of the most convenient and frequently used techniques going beyond 
Poincarc asymptotics is summation to the least term. When x is very large, 
the terms of yo start by decreasing rapidly (klx~ k ~ 1 <C (k — l)lx~ k if k <C x). 
It is tempting to keep adding these terms as long as they continue to decrease. 
The least term, for k = k(x) = [\x\\, evaluates to k(x)\/x k ^ +1 oc x _1 / 2 e _ ' x '. 
Even taking into account the slight ambiguity as to where to stop adding the 
terms, the numerical precision attained at this stage appears to be high enough 
to measure the small exponential beyond all orders. It turns out that in any 
direction towards infinity there exists (see e.g. Costin and Kruskal 1996) an 
actual solution of the equation which lies within 0(x^ 1 ^ 2 c~^) of the optimally 
truncated yo ( a nearest solution); this solution, yo-, is then necessarily unique 
for any fixed arg(cc), since different solutions are separated by Const. e~ x . 

As a manifestation of the Stokes phenomenon, the association between j/o 
and j/o changes with £ = argx. In our example yo — e~ x Ei(x) ±ttic~ x if ±£ > 
and yo = e _x Ei(a;) when £ = 0. 

Since the solutions yo do not have singularities for large x, nothing abrupt 
should happen to the terms beyond all orders, either. Michael Berry discovered 
that there is indeed an intermediate range argx ~ |a;| -1 / 2 where the constant 
C changes smoothly, in a surprisingly universal way. 

We prove that Berry smoothing applies to all nonlinear systems satisfying 
(al) and (a2), and the transition function is the same for all these systems. Tech- 
nically, we show that Berry smoothing is governed by the dominant behavior at 
the singularities in Borcl space nearest to the origin, and these singularities have 
the same essential features for all equations within our class. Going beyond non- 
resonance, the picture changes. However, properly interpreted, Berry smoothing 
might persist in many resonant cases. See Appendix 2 for a discussion. 

As in the toy model above, in some cases (linear homogeneous second order 
and first order equations notably, see Costin and Kruskal 1996) there is only 
one solution of the original differential equation nearest to the divergent series. 
In general there is a restricted family of such solutions. Olver's (1964) example 
illustrates very well that these special solutions might not be the obvious ones. 
At least for this reason it is important to identify a method of recovering them 
from their transseries. 

Hypcrasymptotics as introduced by Berry leads to a substantial improvement 
in precision over optimal truncation but nevertheless ends up with nonzero er- 
rors. Recent very interesting developments of the original ideas of Berry due to 
Olde Daalhuis (1996), suggest that for some linear differential equations hypcr- 
asymptotics may eliminate all errors, but unfortunately at the price of giving up 
the convenience of optimal truncation by adding up vastly more terms beyond 
the least one in an extended family of expansions. 

At present the only technique which at the same time gives full account 
of the complete formal solutions and is widely applicable is generalized Borel 
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summation (see Ecalle 1981 and 1993 , and also Costin 1995 and 1997 for 
rigorous results in the context of l|1.3p). In this paper we show that there is 
a natural connection between Borel summation and optimal truncation, and 
that optimal truncation is compatible to a unique summation procedure, the 
balanced averaging introduced in Costin (1997). 

1.1 Nonresonance 

Let 9 e [0,2n) and A = (X n , X ip ) where jargA^ — B\ £ (— 7r/2,7r/2) (those 
eigenvalues contained in the open half-plane Hg centered along e l9 ). We require 
that for any 9: 

(1) Ajj , . . . , X ip are Z-lincarly independent. 

(2) The complex numbers in the set {Xi — k ■ A £ Hg : k £ N p , i = 1, 
(note: the set is finite) have distinct directions. These are the Stokes lines dj ; k- 

That the set of A which satisfy (1) and (2) has full measure follows from the 
fact that (1) and (2) follow from the condition: 

(m, m' £ Z n , a £ K and (m - am') ■ A = OJ =>■ = can'j (1.5) 

Indeed, if (|1.5fl fails, one of 'SiXj, QXj is a rational function with rational co- 
efficients of the other 5fAj and SAj, corresponding to a zero measure set in 
R 2 ™. 

Notes, (i) Since only small exponentials arc allowed in a proper transseries, 
our conditions only restrict those (sets of) eigenvalues which are associated, 
through (|1.4|) , to exponentials that can be simultaneously small in some direc- 
tion. 

(ii) If condition (1) fails, then higher terms in the transseries (which are 
present for nonlinear systems) may resonate, that is k • A = k' • A for some 
k > 0,k' > (k 7^ k'). This affects the equations characterizing yk- However, 
we think that our results should not change in any substantial way unless in fact 
Xi = Xj for some i ^ j. See Appendix 2 for examples of this latter situation. 

(iii) In the Airy and PI equations discussed in the introduction we have, after 
normalization, Ai,2 = ±1 and the assumptions of nonrcsonnance are satisfied. 

(iv) For linear systems, a weaker nonresonnance condition would suffice for 
our analysis, namely that argAi = argAj for i ^ j (see Costin 1995). This 
condition, as well as (2) above, insures that the Stokes directions are distinct. 

1.2 Normalization and notation. 

By means of normal form calculations (changes of variables), assuming 

(al) and (a2), it is possible to prepare (|1.3|) so that (see Wasow 1968 and 

Tovbis 1992) 

(nl) A = diag(Ai) and 
(n2) B = diag(A) 
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(To this end, it is not necessary that the original matrix B in (|1.3fl is diagonal- 
izable) . 

For convenience, we rescale x and reorder the components of y so that: 
(n3) |Aj| > 1 for i — l...n, Ai = 1, and <f>i < cf>j (cf. (a2)) if i < j, where 
4> = arg(Ai). (To simplify notations, we will formulate some of our results 
relative to Ai; they can be easily adapted to any other Ai, |A^| = 1.) 

To unify the treatment, by taking y = yix~ N for some N > 0, we ensure 
that 

(n4) K(/3 j )<0, j = l,2,...,n. 
(There is an asymmetry at this point: the opposite inequality cannot be achieved, 
in general, as simply and without violating analyticity at infinity.) Finally, 
through a transformation of the form y <-> y — a kX~ k we arrange that 

(n5) f = 0{x~ M - 1 ) and g(x-,y) = 0(y 2 , sc _M_1 y)- We choose M > 1 + 
maxi 5f(— (3i). 

Under the assumptions (al) and (a2) and after the preparation (nl) through 
(n5), the system l|1.3|) admits an n-parameter family of formal exponential 
series solutions, (|1 .4|) (see Iwano 1957 and Wasow 1968; in Costin 1997, a brief 
derivation in the context of Ijl.^Jl is given). In l|1.4(l . C G C" is an arbitrary 
vector of constants, and m.; = 1 — [A J ■ 

The terms of an asymptotic expansion are well-ordered with respect to <C. 
Thus, (|1.4|) is asymptotic in a given direction iff any ascending chain k! • 
Xx) < 3fJ(— k 2 • Xx) < . . ., kj ^ kj, in (|1.4|l is finite (agreeing to omit the terms 
with Ci = 0). For (|1.3fl the condition reads: 

(n6) ^ + <j>i e (-7r/2,7r/2) for all i such that d ^ 0. This selects the 
eigenvalues that lie in a half-plane centered at x. Without loss of generality, 
we let arg(x) vary around arg(Ai) = 0. From now on, A = (A^, . . . , Aj ni ), 
(3 = ([3^ , (3 i2 , . . . , (3i ni ), m = (mjj , m i2 , . . . , rn ini ) and (3' = (3 + m where the 
index is selected by (n6). 

We will henceforth consider that (|1.3|) is presented in prepared form, and 
use the designation transseries only for those formal solutions satisfying (n6). 

With respect to Borel summation we are using the results and notation of 
Costin (1995 and 1997). (A summary is presented in Appendix 1.) 

2 Main results 

All the analysis in this paper is based on the structure of singularities in Borel 
space of the solutions of l|l.l|) and does not otherwise use l|l.l|) . The results and 
proofs can be easily adapted to other functions than solutions of differential 
systems, having a similar singularity structure. 

To simplify the notation, as mentioned in the introduction, we formulate 
all results relative to Ai = 1; they can be reformulated without any difficulty 
relative to any other A; of modulus one. We can reduce the discussion to the 
case when 

Si ^0 
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(cf. theorem [3 see ^2 .21 for a discussion). 

(a) Behavior for large r £ N of the coefficients a rj o in the asymptotic series yo. 
Let 

r(r + 

Proposition 1 For r — > oo wii/i r£N me /icwe 

^1,0 = ^(0)= ]T (e^ + h^r))^ (2.2) 

i;|Aj|=i 

where hj(r) ~ r ~ 1 X^fclo n j; fe7 * _fc ^ or ^ ar 5 e r - ^ e minimum of |a ri0 x _r | zs 
reached for \r — \x\\ < Const, where the constant does not depend on x, r. 

Remark. Wc can write a,-.o = -|=i a r,0;[j] with a^y] ~ ejlYj and 

thus, defining y .yj = J2 r a r,0;[j] x ~ r we nave 

The terms of each series yo;[j] have the same argument (phase) iff a; belongs to 
the j-th Stokes line, i.e. iff arg(x) = <pj. This is a precise formulation of Dingle's 
rule of signs (cf. Dingle 1973, pp.7). For resonant differential systems the rule 
needs to be reinterpreted. See also the Appendix for some further comments. 

(b) Difference between solutions and their optimally truncated asymptotic se- 
ries. 

Theorem 2 (i) Let y be a solution of fl.lty and x = re 1 ^ , with r £ N and 
£ £ [0,27r). Then the difference between y and the optimal truncation of yo is 
of the order of magnitude of the least term, 

y(x) - y [ Q ] (x) = 0(& a ^ r r~ r ) (for fixed £, as r -> oo) (2.3) 

if and only if: (1) y is the balanced Borel sum (see Eq. \4-9^ and theorem\B$ of 
the formal solution y , i.e., 

y = CBiyo + £ cV^aT^Bi y k (2.4) 

|k|>0 

and (2) Cj — for all j such that \Xj\ = 1. 

(wj ^4 similar estimate holds for the higher series in the transseries: For any 
k and large r we have 

CBiy k (x) - y£ ] = 0(a k , r r" r ) (as r -» oo) (2.5) 
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(c) Berry smoothing. We now focus on the thin parabolic region (the Berry 
scale) |ar| » 1, arg(x) — arg(Aj) = fi|a;| _1 / 2 near a Stokes line, say R + (the 
Stokes line for Ai). 

Theorem |3 shows that on the Berry scale near a Stokes line, the constant C 
("beyond all orders") of a given solution changes, and the transition is given, 
as predicted by Berry, by an error function. 

Theorem 3 Let x = re' ar ' , with Q S M. fixed and reN; as r — ► oo we have 

O3iy (x) - yM = ±Si erf (J^J y ei + oie^x^') = 

ift erf (^-^Q e- x a;- /3 'ei + o^-^') (2.6) 

We note here that there always exist directions along which the smoothing 
is visible (namely the directions fa such that |Aj| = 1 in our normalization). For 
other Stokes directions, O(a.o^ r f r ) may be much larger than the contribution 
of erfc. For the directions of Xj with |Aj| > 1 the smallest term of the series is 
too large for the Berry transition to be seen. The change in the terms beyond 
all orders when crossing the Stokes line is always measurable by Borel summa- 
tion (cf. theorem [5] (ii) and Costin 1997) but might be too small for optimal 
truncation. 

(d) Connection with Borel summation. 

In the context of generic nonlinear systems there is a one parameter fam- 
ily of distinct generalized Borel summation formulas which are compatible with 
all usual operations (addition, multiplication, differentiation, composition and 
their inverses, when meaningful) and which preserve reality and asymptotic 
inequalities (see Costin 1997). In other words, the summation operators can- 
not be intrinsically distinguished, in terms of their algebraic and/or classical 
asymptotics properties. 

Nevertheless, proposition^Jshows that only the generalized Borel summation 
with a = 0, the balanced average, is compatible with least term truncation. 
Balanced average is thus the only Borel summation method that guarantees 
optimal least term truncation properties (incidentally, the balanced average is 
also the most symmetric and "natural" in its family). We note that in our 
context it can be shown that the balanced average equals the (apparently) more 
complicated median average of Ecalle. 

Proposition 4 Let y = £B a yo. Then 

y(re i? ) - y l o ] (re*) = 0{^ r r- r ) (2.7) 

for large r and all £ if and only if a — 0. 

Note: with yo fixed, y is a different solution in different Stokes sectors, see 
theorem El 



2.1 Proofs and further results 

To simplify the exposition, detailed proofs are provided for the principal series 
yo- There are only minor adjustments necessary to cover all other components 
yk of the transseries; we mention them at the end. The notation is the same as 
in Costin (1997) and is explained in AO We let Y = Y . 
Taking 

y = jY{p)e- x Up (2.8) 

where C is a contour in 1Z such that $i(xp) > (or a linear combination of such 
contours as in the balanced average), we get by successive integrations by parts 

y(x) = Y W (0)x- fc - 1 + x- r ~ l I Y (r+1 ) (p)e- px dp (2.9) 

Taking x = re~ ltp we get 



y( x ) ~ E ^-ITT- = r-' r - 1 e- i ('-+ 1 ^ f Y^'' +1 ^ (p)e _pre ~*dp (2.10) 

k=0 X J C 

Equation (|2.1U|) represents the error in summation near the least term; we will 
show that for large r 



J Y (r+1) (p)e- pro_i *dp = O (Y (r '(0)) (2.11) 

2.2 Asymptotic behavior of Y^ r ^ for large r G N 

(1) If, exceptionally, all Stokes constants are zero (i.e. Sj = (cf. I|4.5|l 'l for 
j = 1, . . . , n), then Y is entire and |Y(p)| < exp(z/|p|) for some v (theorems[3 
and contour, 



for large r. This implies at once that y IS convergent in l/x. This case is thus 
trivial: summation to the least term is nothing else than convergent summation, 
and the sum is a true solution of the equation. 

(2) By theorem|7|for all i such that Si = 0, Y is analytic along R + Ai. We thus 
only count the Ai with Si ^ and rescale variables and change notations so 
that min{|Ai| : S t ^ 0} = 1 and /3 = 5R(/3i) = mm{3?(A) : |A<| = 1, Ss ^ 0}. 



<lL ex P(^) 2nr~V2^e^-V (2.12) 

27T 7 ,r + i 



9 




CO 
+ 



Fig. 1 

Proposition 5 For large r e N we have 



YM(0) . E /*"~> ^ + 0(r!(1 + „ 13) 



Proof. Let m = 1 + max^m^. Then ^'"Yfjx: 1 ^) is continuous in <f> in 
[0, 27r]\{0i, i = l...n} and has at most jump discontinuities at {4>i,i = l...n} 
(theorem®. Writing 2?ri(r + to)!YM (0) = § V m Y(s)s- r - m - 1 ds, pushing the 
contour of integration towards the circle of radius 1 + e (Fig. 1), and letting 
K = max P m Y((l + e)^) and N x = #{i : Si ^ 0} we have 



Iv m Y(s)s- r - m - 1 ds-Y^ f V rn Y(s)i 

J i=l 

< 2TrK{l + e)- (r+m+1) 
Using theoremElwe have, with E\% = 0(H(1 + e)~ r ), 



l ds 



(2.14) 
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YW(0) = [m . . ' J] / ' SjV™-™^^ (s - A J -)s- r -™- 1 ds + Si 

71-1 j-=i •'Aj 

= E Ci A ^■Y je3 ( S -A J ) S -^- 1 ds + J B 2 (2.15) 



where the j-th term has been integrated by parts m — rrij times and we used 
theorem[7|to show that the boundary terms vanish at s = Xj, and are 0(r!(l + 
e) _r ) at + e). By theorem (ii) we have for small p 



Y jej (p)=pV- 1 A j (p) (2.16) 

and Aj(0) = ej/T(j3'j). Changing variables to s = Xj(l + z) in l|2.15(l we get, by 
Laplace's method, 



2tti j J 

= ^^M^H (2.17) 



2?ri(r + rrij + lpA*-" 



where Gj(z) = O(z) are analytic and O(z) for small z and thus the Hj(r) 



have an asymptotic series in r 1 , and Hj(r) = 0(r 1 ). Combining all these 



estimates we get 



yW W = £ SS^r^ + Ofr- 1 )) (2.18) 



, 27rie i ( r+1 -ft 



so that, with x = re 1V , 



yW(o) 



-r-1 i( 



V S 3 r ( r -0J + 1 ) ( e . , (i/ r )) (2 191 



x r+1 ^ 27rie 

j';|Ail=i 

Because of the e^, the factors inside the sum cannot cancel each-other. The 
rest of proposition [5] is immediate. 



3 Asymptotics of Laplace Integrals 

We distinguish two cases: 
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3.1 Regular directions 



If arg(x) = £ then, by convention, the direction of the Laplace transform integral 
is argp = ip = — £. 

If ip ^ 0j then Y is analytic and exponentially bounded along argp = ip 
(theorems □ and 0. We may assume ip G (0, fc) (cf. fl2~2|) . Let M G N be large, 
subject to the following conditions: 

(ci) ip± arcsin(M _1 ) G (0,<£ 2 )- 

(c 2 ) dist({fc^ : t > M}, {RAi}<=i.. n ) > a > 1 

(c 3 ) M cos(^) > 2. 

Let 5i > 2e be such that M + Si ^ {N|Aj| : i = 1,2, ... ,n} (in particular, 
there are no singular points of Y on the circle of radius M + Si). Let K\ = 
max t < A/+5l \Y(te lip )\ and K 2 = maX| s | =M+(5l |Y(s)|. 

Let N = M + By theorem □ Y(p) is analytic in J M = < 

M + <5i,arg(p) 7^ <pi\ and by theorem (i), j> N \ m \Y is uniformly continu- 
ous in the interior of Jm\ let K 3 = max\ s \ = M+8i;j<\ia\ \P J Y(s)| < (A'i + 
i^ 2 )(2iV|m|) w l m l. By theorems and (ii) we can choose v large enough so 
that for arg(p) ± arcsin(M ) G (0, fa) 

|Y(p)| < exp(KN + 1)) (3.1) 

We note that M, ^, K3, N, Si are chosen independently of r and of x. 
By (c 2 ) and (|3.1|l for t > M we have 

|Y«(*e^)|<^exp(K* + 2)) (3.2) 
For |x| > Mv(M - 2) _1 we have 

/ e- p:i; Y (r) (p)dp 

where K4 depends only on v and M . By construction, for t < M we have 
\te iv - s\ > 2e if |s| = M + <5i. We then have (sec Fig 1.), 



< 



2tt v — \x\ 



<"-W> < K 4 rle~ 2 ^ 



(3.3) 



(r) , _ (jV|m|+r)! t V^Y(s) 
1 } 2ni ^(s-ufry+NM+i 



(N\m\ +r)\ 
2~7ri 



E 



pJV|m|Y( s ) 



ds + r\e~ 2r E 2 



c . (s - tc i f) r+Nm '+ 



T ds + rlc- 2r E 3 (3.4) 
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where 



fi = (Nmi + r)\/(2Tri) 



Ei contains the contribution of the outer circle, and in the last equality we 
integrated by parts |m| - m t times and included the boundary terms in £3. 
We have \Ei\ + \E^\ < E4 for some E4 depending only on K$,M and |m|. In 
conclusion, 



/ e- px Y ( - r) (p)dp = Y j f J e~ px dp 
Jo - =1 Jo Jc 



Y(s) 



ds 



+ r\(e- 2r + e~ 2lxl )E 5 =J2fj / dsT m i N Y{s) 
3=1 Jc i 



Me"* 



e px dp 



(s - pe l v) Nm J+ r + l 
+ rl{e~ 2r + e- 2 ^)E 5 (3.5) 



for large enough r, a; where E5 is independent of r, x. 
We have, with p = tc iip , x = rc~ iip , R = mjN + 1: 



Mc" 



- px 



(s - p) 



dp 



- px 



The behavior for large r of the last integral is obtained by standard steepest 
descent: The function f(t) = c _ *(k — t) -1 has a saddle point at t = k — 1 and, 
for k ^ [0, 00] the steepest descent path originating at zero starts in a direction 
opposite to the pole k and continues to +00. For arg(,s) 7^ ip, the contour in 
l|3.6J> can thus be deformed to a steepest descent curve without crossing the pole 
or the saddle. The main contribution to the integral comes then from a region 
where t is small. We have the estimate for (|3.6H : 



o (s - te^y+R 
1 



dt 



cxp (-(r + R) ln(s - tc iip ) - tr)) dt 



S r+R 



exp(-t(r - (r + R)e lip s~ 1 ) + 0{r 2 ))dt 

= [rs r+R (c iv s^ 1 - l)]"^! + Oir- 1 )) (3.7) 



As p Nm iY ± are continuous along Cj, we obtain 



(M+<5i 



\j(l+e) 



„-r-R • 



>Nr, 



'Y±(s) 



e i V g-l _ 1 



ds 



<(l + e) 



-r-R 



(3- 
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for some Eq independent of r and x. Since r p m iY ± are continuous for small 
e on the segment between Xj and Xj(l + e) we have, using (|4.13|) . (|3.8|l and 
integration by parts, with Cj being the part of Cj within e distance of \j, 



f 3 f dsP m * N Y(s) f 
Jd Jo 



Me " e- px dp 



l (s-pe i »') Jv ^+ r + 1 
= fj / dsV m > N Y(s) / . .„ P - + 0((1 + e)" r ) 

= / / pmiY^C-P)" 1 ^^- 1 ) / . . P - +Q((l + e)- r ) 

A 3 (l+e) fMe* „-px Hn 

S.-Ye.fs-A,-) / t ■ , t xi +0((l + e)- r ) 

JJ J (s -pe 1 ^) m ^+ r + 1 u ; ; 

ds t J °V ■ f^— (l + 0(r- 1 )) + 0((l + e)- r ) (3.9) 

with 

e, = (m, +r)!/(2iri) 

The last integral is very similar to (|2.15|) and is estimated in the same way, 
giving 



.V,Y,. (s A, 



d 



rs r+m J+ 1 (e i ^s- 1 - 1) 



where, in view of l|2.1(J|) we need to take r + 1 instead of r. Noting that (r + 
rrijy.r-ft- 1 = T(r - (3j + 1)(1 + O^ 1 )) and comparing (|3~TU|) to QJg) , (j!TTTJ) 
is proven. 

□ 



3.2 Stokes directions; balanced averages; Berry smooth- 
ing 

This section deals with the special and important case of exponential asymp- 
totics on and near the Stokes line corresponding to one of the Xi of largest 
module (|Aj| = 1 in our normalization; without loss of generality we take i = 1). 
Berry's smoothing formula is proved, and we show that the balanced average 
is the only Borcl summation process compatible with optimal truncation. The 
generalized Borel summations that have good algebraic and analytic properties 
are given by l|4.9|l . Using theorem [7| (i) which shows that for p G (0, 1) we have 
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Yq~(p) = Y (p), formula l|4.13|) . and the estimates in theorem^ (i) for (£Yk)~ 
we see that 



y := y a = CB a y Q = (1 - a)CY+ + aCY^ + O^ 2 -^) (3.11) 

Choosing < 5 < mini|arg(Ai) — arg(Ai)|, relation (|2.1U|> reads, in view of 
theorem [H] 



y Y«(0) 

where f. . means 

•j <Q:a> 



e ^ x ) 



7 n J <a;oo> 

(3.12) 



ac 



a / +(1 — a) 



In order to prove theorem and proposition 01 we need to estimate Ij3.12|l near 
the Stokes line. Replacing everywhere by J <a . M> and taking y> = 0, the 

calculations leading to l|3. 5J) work without any other change and we get: 



f f f •p TO 3 

/ e-^Y^(p)dp = y2f j e~**dp -—- 

J<a;oo> =1 J<a:M> JCj [ S ~ P) 

+ r\E 5 (e- 2r + e-W) = fj [ dsV m > N Y(s) f 



Nrrij+r+1 



ds 



e px dp 



( S - p) Nm i+ r+1 



>Cj J<a-M> 

+ r\E 5 (e- 2r + e~ 2 ^) (3.13) 
For j ^ 1, s — p does not vanish for s £ Cj, | arg(p)| < (5, and thus 



f> / dsr^ N Y( s ) [ 

3=2 J Cj J <' 



e~ px dp 



ICj J<a;M> ( s V) 1 

and the estimates leading to (|3.1l)|) are valid if we take tp = Q. We get, for j ^ 1, 



M e- pa: dp 
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and 



^(l+e) S i Y e .(s-X i ) 



s,(r :^,„- 1 ( 3 .i6) 



2mr{r + mj + 1) J \ J ] (Xj L - 1) 
We are therefore left with the problem of estimating 

e~ px dp 



J:= / dsT miN Y{s) 

JCt J<a;M> 



(s — p\Nmi+r+l 



(3.17) 



We consider the case f2 < 0, the other case being treated symmetrically. 

Proposition 6 Let 6 G (0,tt/2), n G N, a; = re i0r ~ 1/2 , < 0, s > 1 and 
R = r + n. 

i) For large r G N 



9a = a 



r -l/2 



Me" if 



e p:r dp 
(s — 



-^ 2 + 



f2 r — icr H — (erfl r + ir^ 

27ri 



n dt + 0(r" 



-(!-«)- 



O r = -iVF(e i0 ^ 1/2 
iij When a is small we have 



(R-iy. 

1); tr = VF(«-l) 



(3.18) 



(3.19) 



g a = ie 



( erf ( — + 1 ) E 1 - 2(1 - a)E 2 



-n 2 /2 0(r -i 



(3.20) 



where E h2 = (1 + O^- 1 / 2 , ct)), erf (x) = 27T- 1 / 2 /* exp(-t 2 )di. 
Proof. 



Mc e -px dp X r-o e ~psx dp 



(s- P y 



{1-P) H 

1 / f 00015 e- psx dp 
^ [Jo (l^F 



+ 0(e~^'°)\ (3.21) 
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Furthermore 



c-P x d P i Jrl(lqFl) h-x x ._ /"°° e~ pa dp _, f 00 e-^dp 



7o 7o (s-ip) R 

(3.22) 

Taking a = 1 /2 — e for small e we have 
and 

/" e -isxt-Rln(l-it)^ _ C e -isxt+Rit~Rt 2 /2^ t + (9(1) 
JO Vo r 

e _ isx4+ffit _ flt 2 /2dt+0( ^ = r _ 1/2 jr-V*z(R-s X )-£* 2 dz + (i^ 



(3.24) 

Part (i) of proposition El follows by combining (|3.19l . l3~2"4*|l ). Part (ii) is a 
straightforward calculation from (i), using Stirling's formula and Laplace method. 

□ 

We note that with f2, 7*1 held constant the integral in H3.18|) is bounded for s > 1 
uniformly in r G N. Thus, for s > 1 + e we have g a < ii^e - " for some -K2 and 
thus we obtain from (|3.17() . proceeding as for (|3.9() . 



d S P miJV Y( s ) / 6 - + £ 3 



Ci;|s|<l+e J <a-M> 



(3.25) 



where £^4 = 0(e re ), and /i(s) is smooth on [1, 1 + e]. Taking <5 small and r 
correspondingly large we have, in view of proposition (i), 



17 



fl Jl J< a ;M> {s-p) mi+,+1 

(3.26) 

where E 5 = 0(c~ r ). Using now proposition[S](ii) to evaluate the inner integral 
we obtain 



+(erf(n/-\/2) + 1)^6 - 2(1 - a)£ 7 + e _na/2 0(r _1 )) (3.27) 



where £^,7 = 0(r -1 / 2 ). where, as in (|3.1U|) we need to take r + 1 instead of r 
(cf. i|2.11[) '). Taking A = 1 in the expression l|3.10|l . which was shown to be of 
the order of the least term, we see that the least term is 0(r" 1 / 2 ) smaller than 
[ET27|l unless 

erf(fl/V2) - 1 + 2a = (3.28) 

(a) Asymptotics along straight lines: If SI = (arg(x) = 0) we see that only 

a = 1/2 

(corresponding to the balanced average, which on the interval (0, 2) is the half 
sum of the upper and lower continuations) ensures errors in optimal truncation 
of the order of the least term. 

(b) Asymptotics along parabolas and Berry smoothing. Theorem[2]is straight- 
forward application of Stirling's formula to <|3.27l) . (Note also that by theorem^ 
(ii) and theorem (ii) y ei ~ e~ x ei.) Convergence of the Puiseux series 
near the singularities of Yo in Borel space was the key in proving universality 
of the Berry transition, since using this convergence, the calculation reduced in 
effect to the case of one particular function, Ei(cc). 

The same line of proof as for yo works for yt as well. Indeed, we note that 
have the same singularity structure in Borel space as yo- This fol- 
lows from theorem (ii) and the fact that convolution with a locally analytic 
function preserves the convergence of Puiseux series (cf. Example (3a) in ap- 
pendix 1). □ 
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4 Appendix 1. 



4.1 Results on Borel summation 

In this paper we make use of some of the results in Costin (1997); for convenience 

we summarize them below. 

We use the convention N = N U {0}. Let 

W = {pe C :p^ k\ t ,Vfc G N,i = 1,2,.. .,n} (4.1) 

The directions dj = {p : arg(p) = 4>j},j = 1,2, ...,n are the Stokes lines 
(Note: sometimes known as anti-Stokes lines!). We construct over W a surface 
7?., consisting of homotopy classes of smooth curves in W starting at the origin, 
moving away from it, and crossing at most one Stokes line, at most once (Fig. 
2): 

U := {7 : (0, 1) W : 7(0+) = 0; ^It(*)I > 0; arg( 7 (i)) monotonicj 

(4.2) 

The Laplace transform along a direction <^> of a function F C^F will depend 
in general on </>; the usual convention is to choose <f> so that xp G R + . Thus, 
the Borel sum of / in the direction x, if it exists, is defined as CM x \Bf with 
4>{x) = — arg(x). By (n6) and the agreed association between p and x, and 
since Laplace integrals will not depend on the direction of p until one of the 
Stokes lines is crossed, we may as well assume that the direction of integration 
is either dj or is arbitrarily close to it. Define TZ\ as the restriction of 1Z to 
arg(7) G (^„ — 27r, 1P2) where ip n = max{— w/2, 4> n — 2ir} and ^2 = min{7r/2, ^2}- 







^ 












































• 3 ». 









Fig 2. TTie paths near A2 belong to 1Z. 
The paths near Ai relate to the balanced average 
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We denote the analytic continuation of / along a curve 7 by AC 7 (f). For the 
analytic continuations near a Stokes line di we use notations similar to Ecallc's: 
/~ is the branch of / along a path 7 with arg(7) < <f>i, while f~ k+ denotes the 
branch along a path that crosses the Stokes line between kXi and (k + l)Xi- We 
use the notations V f for f Q f(s)ds and V 7 f if integration is along a curve 7. 

We write k > k' if h > k[ for all i and k y k' if k > k' and k^k'. The 
relation >- is a well ordering on N ni . We let e_,- be the unit vector in the j th 
direction. 

Formal expansions are denoted with a tilde and capital letters Y, V . . . will 
usually denote Borel transforms or functions otherwise associated to Borel space. 
For notational convenience, we will not however distinguish between Y& = By^, 
which turn out to be convergent series, and the sums of these series Y k as germs 
of ramified analytic functions. 

We have 

s(x,y) = gi^y 1 = J2 gsa^y (M > x , |y| < yo) (4.3) 

|1|>1 s>0;|1|>1 

where y 1 = y^ 1 • ■ • y l ™ and |1| = l\ + ■ ■ ■ + l n . By construction g s .i = if |1| = 1 
and s < M. 

Theorem 7 (i) Yo = By® is analytic inlZU {0}. The singularities of Yq 
(which are contained in the set {IXj : I G N + , j = 1,2, ... ,n}) are described as 
follows. For I € N + and small z 



Y±(z + lX J )=± 



(±S J ) l Hz)°' 1 Y le] (z) 



(hrij 



+ B w (z) (/ = 1,2,...) (4.4) 



where the power o/ln(z) is one iffiPj G Z 7 and Ajj, Bjj are analytic for small z. 
The functions Yk are, in addition, analytic atp = IX j, I £ N + , iff, exceptionally, 



S j =r j T(p' j )(A hj ) j (0) = 



(4.5) 



-27ri otherwise. The Sj are Stokes 



where rj = 1 - e 27 ^" 1 ) iflfy <£ Z and rj 
constants, see theorem MCX 

(ii) Yk = Syk, |k| > 1, are analytic in 1Z\{— k' • A + Xi : k' < k, 1 < i < n}. 
For I G N and p near IX j, j = 1, 2, ... ,n there exist A = Akj; and B = Bkj/ 
analytic at zero so that (z is as above) 



Y^z + lX,) 



± 



ukj + l 
JJ [ l 



^'+^- 1 (lnz)°' 1 A k;j (.) 



ln(z)°- 1 Y k+/ej (z) 



(Irrij ) 



ZB k y (z) — 

m^{z) {1 = 0,1,2,...) (4.6) 
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where the power of hi z is iff I = or k • (3 + 1/3 j ■ — 1 ^ Z and A^j = e j /^Wj) ■ 
Near p € {— k' • A : -< k' < k}, (where Yq is analytic) Yk, k / have 
convergent Puiseux series. 

Let Z?yk be extended along dj by the "balanced average" of analytic contin- 
uations 

By k = Yt = Y+ + £ - (Y k - Y k (4.7) 

i=i 

The sum above coincides with the one in which + is exchanged with — , account- 
ing for the reality-preserving property. Clearly, if Yk is analytic along dj , then 
the terms in the infinite sum vanish and Y k Q = Yk; we also let Y k a = Yk if 
d ^ dj, where again Y^ is analytic. It follows from i|4.7[) and theorem [3] below 
that the Laplace integral of Y k a along R + can deformed into contours as those 
depicted in Fig. 2, with weight 2~ k for a contour turning around (k + 1)X\. 
More generally, we consider the averages 

oo 

S Q y k = Y£ = Y+ + £ cJ (Y k - Y k (4.8) 
i=i 

and correspondingly 

(LB) a y k := CYZ (4.9) 

With aeR, this represents the most general family of averages of Borel sum- 
mation formulas which commute with complex conjugation, with the algebraic 
and analytic operations and have good continuity properties (see Costin 1995). 
The value a = 1/2 is special in that it is the only one compatible with optimal 
truncation. 

Theorem 8 (i) The branches of (Yk) 7 in 1Z\ have limits in a C* -algebra of 
distributions, T>' m y (K + ) C T>' (cf. $ \4-ty Their Laplace transforms in T>' m „(R + ) 
£(Yk) 7 exist simultaneously and with x € S x and for any S > there is a con- 
stant K and an x\ large enough, so that for ^t{x) > X\ we have |£(Yk) 7 [x)\ < 

In addition, Y^pc 1 ^) are continuous in <fi with respect to the T>' m v topology, 
(separately) on [ip n — 2tt,0] and [0,ip2\- 

If m > maxj;(?7z,;) and I < mini|Ai| then Yo(pe 1<ji ) is continuous in <f> £ 
[0, 27r]\{(^j : i < n} in the T>' m „(R + , I) topology and has (at most) jump discon- 
tinuities for <f> — (f>i. For each k, |k| > 1 and any K there is an I > and an m 
such thatYkipe 1 ^) are continuous in (j) £ [0, 27r]\{</>.;; — k'-A + A; : i < n, k' < k} 
in the T>' m „((0, K), I) topology and have (at most) jump discontinuities on the 
boundary. 
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(ii) The sum \4-T\) converges inT>' m v (and coincides with the analytic contin- 
uation o/Yk when Yk is analytic along M. + ). For any 8 there is a large enough 
x\ independent of k so that Y£ a (p) with p G 1Z\ are Laplace transformable for 
$t(xp) > X\ and furthermore \{CY^){x)\ < (S' k '. In addition, if d ^ R + , then 
for large v, Yk G L\{d). 

The functions £Y£ a are analytic for ^ft(xp) > x\. For any C G C™ 1 there is 
an xi(C) large enough so that the sum 

y = CY b Q a + C k e- k Xx x- k f3 £Yl a (4.10) 

|k|>0 

converges uniformly for ^R(xp) > xi(C), and y is a solution of \l.£fy . When the 
direction of p is not the real axis then, by definition, Y^ a = Yk, C is the usual 
Laplace transform and \4-lU\j becomes 

y = £Y + C k c- k Xx x- kl3 CY k (4.11) 

|k|>0 

In addition, £Y^ a ~ yk for large x in the half plane di(xp) > x\, for all k, 
uniformly. 

Hi) The general solution of 11.31) that is asymptotic to yo for large x along 
a ray in S x can be equivalently written in the form \4-lty or as 

y = £Y± + J2 C k c- k A2; a;- k/3 /:Y± (4.12) 

|k|>0 

for some C (depending on the solution and chosen form). With the convention 
binding the directions of x and p and the representation form being fixed, ( cf. 
the beginning of j O . l\j ) the representation of a solution is unique. 

Theorem 9 i) For all k and 5i(p) > j, 3(p) > as well as in T>' m v we have 

*F T <P) Y^-^ip) = (±5iy / j ) (Y£ +jei ( P j)) ^ (4.13) 
and also, 




(4.14) 



ii) Local Stokes transition. 
Consider the expression of a fixed solution y of \1.$) as a Borel summed trans- 
series \4.1U\ ). As arg(x) varies, \4-Mty changes only through C, and that change 
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occurs when the Stokes lines are crossed. We have, in the neighborhood of R + , 
with Si defined in j4-5\/ -' 

( C- = C(-0) forZ<0 
C(0 = { C° = C(-0) + |Siei /or£ = (4.15) 

[ C+ =C(-0) + 5iei /or£>0 

Remark 1 In view of \4-l£fy the different analytic continuations of Yo along 
paths crossing K. + at most once can be expressed in terms of Y jei . The most 
general formal solution of \ l.ty) that can be formed in terms of Yj ej . with j > 
is \l-4\l with C\ — a arbitrary and Cj = for j ^ 1. Any true solution of 
based on such a transseries is given in \4-lty with C as above. Any average *4Yo 
along paths going forward in R + such that CAYq is thus of the form 

Theorem 10 Assume only Ai lies in the right half plane. Let 7 ± be two paths 
in the right half plane, near the positive/ negative imaginary axis such that 
\x~P 1+1 c~ xXl | — *• 1 as x — > oo along 7 ± . Consider the solution y of fl.Sp given 
in 14 -lOH with C = Ce\ and where the path of integration is p € M + . Then 

y = (C± I^e^+V^ 1 (i + (i)) (4.16) 

for large x along 7 ± ; where Si is the same as in {4-5$ , [4.15\j . 

Proposition 11 i) Let y\ and y2 be solutions of ll.ty so that yi,2 ~ yo 
for large x in an open sector S (or in some direction d); then yi — y2 = 
Y) j Cje~ Xz i I x~^ l i (e^ +o(l)) for some constants Cj, where the indices run over 
the eigenvalues A^. with the property 5R(A^x) > in S (or d). If yi — yi = 
o(e~ x '3 x ' x~ i ) for all j , then yi = y2- 

ii) Let yi and yi be solutions of and assume that yi — y2 has differ- 

entiable asymptotics of the form Kaexp(— ax)x b {l + o(l)) with di(ax) > and 
for large x. Then a = X{ for some i. 

Hi) Let Uk £ Zr.j for all k, |k| > 1. Assume in addition that for large v 
there is a function S(v) vanishing as v — > oo such that 

su P( T |k| / lUk^e-^l d|p| < K < oo (4.17) 

k Jd 

Then, if y\,yi are solutions of in S where in addition 

yi - y 2 = e ~ X kXxm k I u k (p) cxp(~x P )dp (4.18) 
ik|>i Jd 

where \,x are as in (n6), then yi = y2, and Uk = for all k, |k| > 1. 
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4.2 Focusing spaces and algebras 

The proofs of the properties stated in this section are given in Costin (1997). 

We say that a family of norms HH^ depending on a parameter v E M + is 
focusing if for any / with \\f\\ vo < oo 

11/11, 1 as v T oo (4.19) 

Let £ be a linear space and {||||i/} a family of norms satisfying (|4.19() . For 
each v we define a Banach space £>„ as the completion of {/ € £ : ||/||„ < oo}. 
Enlarging £ if needed, we may assume that B v C £ . For a < (3, H4.19f) shows 
that the identity is an embedding of B a in Bp. Let T C £ be the projective 
limit of the B v . That is to say 

T := |J B v (4.20) 

is endowed with the topology in which a sequence is convergent if it converges 
in some B v . We call T a focusing space. 

Consider now the case when (Z?„, +, *, || ||„) are commutative Banach alge- 
bras. Then T inherits a structure of a commutative algebra, in which * ("con- 
volution") is continuous. We say that (T 1 *, ||||„) is a focusing algebra. 

4.3 Examples 

For more details see Costin (1997). Let K £ M + and iS = SK, ai ,a 2 = {p '■ 
arg(p) G [ai, at2\ C (— 7t/2, 7t/2), \p\ < K} (or a finite union of such sectors) and 
V be a small neighborhood of the origin. V will be the closure of V, cut along 
the negative axis, and together with these upper and lower cuts. 

(1) . L^/C). Let K, = SK,<ji,<j)- The space L 1 ^(/C) with the convolution 
/ * g := p i— * Jq f(s)g(p — s)ds is a commutative Banach algebra under each of 
the (equivalent) norms 

||/!|, = [ e-^|/(*exp(i^))|d* (4.21) 
Jo 

(2) If A' = oo in example (1), then the norms (|4.21|) are not equivalent 
anymore for different v, but convolution is still continuous in l|4.21|l and the 
projective limit of the L^K+e 1 *), ^(R+c^) C i^R+e 1 *), is a focusing alge- 
bra. 

(3a) Tp(S U V). For 3?(/3) > and </>i ^ </> 2 , this space is given by {/ : 
f{p) = P F{p)}i where F is analytic in the interior of S U V and continuous in 
its closure. We take the family of (equivalent) norms 

\\f\\u,p = K sup |e-*/(p)| (4.22) 
sesuv 
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It is clear that convergence of / in \\\\ v ,i3 implies uniform convergence of F on 
compact sets in S U V (for p near zero, this follows from Cauchy's formula). Tp 
are thus Banach spaces and focusing spaces by Q4.22p . The spaces {Tp}p are 
isomorphic to each-other. The application 

(•*■): T Pl x Tp 2 i * T A+ft+1 (4.23) 

is continuous: 

A natural generalization of 7g is obtained taking pi, ... , (3m € C with pos- 
itive real parts, no two of them differing by an integer. If fp = P^-^iip) 
with Ai analytic, then fp = iff A4 = for all i (e.g., by a Puiseux series argu- 
ment). It is then natural to identify the space T{p u ...,p k \ of functions of the form 
fp with ®^_ 1 Tp i . Convolution with analytic functions is defined on T\p ly p k y 
while convolution of two functions in T{p lt p k y takes values in T,p + ^ mod i}' 
We write Tr.x when the concrete values of pi, . . . , pk do not matter. 

(3b) A particular case of the preceding example is A z ,i(S U V) consisting of 
analytic functions in the interior of SUV, continuous on its closure, and vanishing 
at the origin together with the first I derivatives. A z ,i can be identified with %. 

(4) T>' mv , the "staircase distributions". Let T>(0,x) be the test functions 
on (0, x) and T> = 2?(0, oo). Let T>' m C V be the distributions / for which 
/ = Fl\ km) onX>(0,/c+l) with F k £ L^O.fc + l]. There is a uniquely associated 
staircase decomposition, a sequence {Ai(/)} igN = {Ai} igN such that A* £ 
A, = AiXui+i] and 



/ = E A * (mi) ( 4 - 24 ) 

i=Q 

With respect to the norm 

oo 

\\f\\ v , m :=V2j2 vim \\^i\U ( 4 -25) 
»=o 

where |jA||„ is computed from (|4.21|) with K = oo and with convolution defined 

as 



A fe (/ * /) = J2 A * * A ^ - v " 1 { E ( A * * A *Wi] > ( 4 - 26 ) 

(2?' n ,+,*) is a commutative Banach algebra. With respect to the family of 
norms ||||m,i/, the projective limit of the T>' rn v , T m is a focusing algebra. 

For any / £ L^ o (M + ) there is a constant C(y, vq) such that / £ V m v for all 
v > 1/q and 
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Il/lk„„ <C{y ,u)\\f\\ Lla (4.27) 

and formula 1)4. 26[) is equivalent to the usual convolution in this case. 

For a e K + , T>' mv {a,oc) = {/ £ T>' m v : Ai(x) = for x < a} is a closed 
ideal in V' mv (isomorphic to the restriction V' m v (a, oo) of T>' m v to T>(a, oo)). 
The restrictions 2?^ „(a, b) of „ to T>(a, b) are for < a < b < oo Banach 
spaces with respect to the norm (|4.25[) restricted to (a, b). 

The functions in T> (IR + \N) are dense in T>' mv , with respect to the norm 
(H2SJ). 

If we choose a different interval length I > instead of I = 1 in the partition 
associated to 1)4. 24f) . we then write X>£„ „(Z). Obviously, dilation gives a natural 
isomorphism between these structures. If d = {tc 1 ^ : t e M+} is any ray, T>' mv (d) 
and T m ;4> are defined in an analogous way and have the same properties as their 
real counterpart. 

Laplace transforms are naturally defined in T>' m v . 

Lemma 12 Laplace transform extends continuously fromT>(JS, + \N) toT>' m „(K + ) 
by the formula 

°° />oo 

(Cf)(x) :=Y. xmk / e~ SX Ms)ds (4.28) 
fc=o Jo 

In particular, with f,g,h' € T>' m v we have 



C(f*g)=C(f)C(g) 

C(h') = xC(h) - h(0) 

C(pf) = -(£(/))' (4.29) 

For x € SV = {.t : 5R(a;) > f } £/ie sum converges absolutely. Laplace 

transform is, for fixed a continuous functional (of norm less than one) 

on V m v . 

(Cf)(x) is analytic in S u . 

Furthermore, C is injective in T>' m v . 



4.4 Appendix 2. Two examples of resonant equations 

Once resonant equations are allowed, Berry transitions tend to become more 
complicated: y = e~~ x Ei(x) + exp(— x 2 — ice 4 ) is a (least term summable for 
x — > +oo) solution of a linear differential equation with rational coefficients: 



d 
dx 



EL 



2x + 4ix c 



d 2 . ,.d 
X ^ 2+{X + I) d-X 



1 



y = 
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where P(x) = -16.t 7 + 16ix 5 - 4i.T 4 + (4 - 16i)x 3 - 2x 2 - Ax + 1. The term 
cxp(— x 2 — ix 4 ) is not seen by least term truncation on the real axis but becomes 
much larger than the error function contribution before entering the Berry region 
arg(x) ~ |a;| -1 / 2 . 

Borel summation using Ecalle acceleration (see Ecalle 1993) gives an unam- 
biguous description of the solutions of this type of equations and can be used 
to decompose solutions conveniently before analyzing transitions. On the other 
hand, imposing that all Stokes constants are nonzero suffices to exclude this and 
similar examples, but looks rather restrictive and hard to verify. 

We next consider the Berry transition of a family of resonant equations with 
nonzero Stokes constants. The formal solutions of such equations depart from 
Dingle's rule, and also exhibit an interesting splitting of the Stokes rays, with 
two Berry transitions in the same region. (We are mainly aiming at illustration 
and the calculation is heuristic but a rigorous treatment along these same lines 
is not difficult). Let first m £ R and 



(A) Formal solutions; behavior of coefficients. 

oo 

k=0 ' 



Taking y = J2T=o akX k we S et for Uk 



a k+1 = (2k - m 2 )a k - k(k - l)a k -i ( 4 -31) 
With a k = (k — l)\b k we have 



b k+ i = (2 - — )b k - 6 fe _! (4.32) 

which we analyze for k 3> 1 by WKB (an explicit solution is also possible in 
this case). With b k = e Wk we get: 



'/ Til i 1 // 

~ (2 )e w - e w ~ w + * w 

k 

ei»" coshK) ~ 1 - ^(1 + !«/')(! + \w' 2 ) ~ 1 - ^ 



* - T + w " ~ ±i V X T ^ ~ ±l v=k + 5* (433) 



so that 



b k ~ fc 1/4 (yl + c 2imvT + A_c- 2imvT ) and a k = a+ + 

with a± ~ (k - l)L4±ife 1 / 4 e ±2WS (4.34) 
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and, in analogy with the nonresonant case we write y = y + + y~ where ^ = 
SfeLo a k x ~ k - We first see that there is no curve arg(x) = f{\x\) for which the 
terms of (or the terms of any linear combination of y^) have the same phase. 
There are instead two parabolic curves, ±arg(x) = 2m\x\~ 1 / 2 along which for 
n ~ |a;| we have {a^ +1 / x n+1 ] / {a+ / x n ] = 1 + o(l). This does not amount to the 
terms being in phase, but rather is the discrete equivalent of a stationary point. 

To reinterpret Dingle's rule for this example, wc will sec that there exist two 
Stokes parabolas, each associated with one degree of freedom in the original 
equation along which the transitions in the constants beyond all orders, as 
measured by optimal truncation, are maximal. 

(B) Asymptotic solutions of the homogeneous equation. 
Taking y = e w we obtain: 



9 TTi^ / 777,^ 

w" + w' + 2w' + 1 + — = => w' = -1 ± i\ h w" 



w~-x± 2imy/x~ +-\n x=>y ± ~ x V* e -x±X™^ ^ ^ 



as before. 

(C) Berry smoothing. 



Let y = J2k=i a k% k + E± c± ( x )y± ( x )- Thcn Hv] - x 1 S ives 



± 



CI + 2C'+ ( 1 + ?± 

v± 



n(n — l)a„_i (n — l)(n — 2)a ra _ 



.r 



n+l 



x' 



m 2 - 2n - 2 n(n - l)a„_ 1 a T{ 

a n -\ = L + C + + L-C- H 

x n x n ~^~ X r ' 



L+C+ + L-C- 



nl , ( 1 _ b n -i 
b n x 



(4.36) 



with obvious notations. Changing variables to x = nc ll3n above we get, to 
leading order, 



8 l/4 y e -n+iVi(±2m-(!) + ^ 2 Tm3 ( ^ 2mC '± 

^ \ n n 



E 



n 



-iVn/3 



V* A± e ±2imn ' n-V\±m + p) 



(4.37) 



Equating the coefficients of e 2mm 1 we get, with B± = Const + o(l) the system 



2im 



C'l ± ^-C' ± = \B ± n-^-^ ±m ^{± m + /?) 



(4.38) 
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where we change variables to x = nexp(i/3?i x l 2 ) and get to leading order: 

-C'l(p)±2mC ± (p) =iB ± e-^ 2±mfi (±m + P) (4.39) 
with the bounded solutions: 

C± 09) = i^^erf B± + Const.± (4.40) 
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